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Motivated by second order asymptotic results, we characterize the convergence in law of double 
integrals, with respect to Poisson random measures, toward a standard Gaussian distribution. 
Our conditions are expressed in terms of contractions of the kernels. To prove our main results, 
we use the theory of stable convergence of generalized stochastic integrals developed by Peccati 
and Taqqu. One of the advantages of our approach is that the conditions are expressed directly 
in terms of the kernel appearing in the multiple integral and do not make any explicit use of 
asymptotic dependence properties such as mixing. We illustrate our techniques by an application 
involving linear and quadratic functionals of generalized Ornstein-Uhlenbeck processes, as well 
as examples concerning random hazard rates. 
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1. Introduction 

Let TV {N{B) : B E Z} be a centered Poisson measure over some Borel space {Z,Z) 
and let {gn-n > 1} and {/„ :n > 1} be, respectively, a sequence of real- valued functions 
over Z and a sequence of real- valued symmetric functions over Z^, vanishing on the 
diagonal set {{a,b) G : a = b} . The aim of this paper is to characterize the convergence 
in distribution, toward a bivariate Gaussian law, of sequences of random vectors of the 
type 

{hign)j2ifn):n>l}, (1) 

where 

/ gn{a)N{da), n>l, (2) 
Jz 
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are single Wiener-Ito stochastic integrals with respect to A^and 

Wn)= f f fnia,b)Nida)N{db), n>l, (3) 
J z J z 

is a double Wiener-Ito integral with respect to N (see [24], as well as the discussion 
contained in Section 2, for precise definitions). One of the main motivations for the study 
of random objects such as (2) and (3) is that the (Hilbcrt) space of random variables Z 
having the representation 

Z = c + h{g)+h{f), (4) 

for some real constant c and some (uniquely defined) deterministic kernels g and /, coin- 
cides with the i^-spacc generated by random variables of the type P{N{Ai), . . . , N{Am)), 
m > 1, where P(-) is a polynomial of degree < 2 in to variables (see, e.g., Ogura [14], 
or Rota and Wallstrom [21], for a proof of this fact and for a representation of multiple 
Poisson integrals in terms of Charlicr polynomials). It follows that a characterization of 
the asymptotic normality of a sequence such as (1) is crucial in the study of the asymp- 
totic behavior of random variables that are linear or quadratic transformations of the 
measure N. For instance, we will see below that one can effectively study linear and 
quadratic functionals of stochastic processes constructed from N by first resorting to 
their representation in the form (4). 

The main point in the characterization of the asymptotic normality of (1) is to obtain 
criteria for the weak convergence (toward a Gaussian law) of the sequence {l2{fn) '■n> 1}. 
To do this, we shall make extensive use of the results, involving 'generalized adapted 
stochastic integrals,' developed in [17] and [19]. These results arc based on a decoupling 
technique, known as the principle of conditioning (see, e.g., Xue [27]). In particular, the 
sufficient conditions for the asymptotic normality of {l2{fn)'-n> 1} will be expressed 
in terms of the contraction kernels associated with the functions /„. These contraction 
kernels, whose definition is given in Section 2 below, can be easily computed from the 
analytic expression for /„. 

We will also show that if the variances of hifn) converge to 1 and if some integrability 
assumptions are satisfied (e.g., if the sequence {l2{fn)'^ ■n> 1} is uniformly intcgrable), 
our conditions are necessary and sufficient, and also that they are equivalent to the 
convergence of the fourth moments E[/2(/n)^] toward 3. Note that the number 3 is simply 
the fourth moment of a standard Gaussian random variable: in our opinion, this fact is 
rather striking, as it shows that central limit theorems (CLTs) involving the integrals 
^2(/n), n>l, are essentially determined by the convergence of quantities associated 
with their first two even moments. This implication should be regarded as a drastic 
simplification of the method of moments and cumulants, customarily adopted to prove 
CLTs for polynomial forms in random variables with finite moments of any order. See 
the two surveys by Surgailis, [24] and [25], for more details on this point. 

Our findings extend to the framework of Poisson measures results previously estab- 
lished in [13, 18] and [20] for sequences of multiple Wiener-Ito integrals (of arbitrary 
order) with respect to general Gaussian processes. For instance, in the Gaussian case, 
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the necessary and sufficient conditions for a CLT involving multiple integrals also require 
the mere convergence of second and fourth moments and can cquivalently be expressed 
in terms of contraction kernels. Note that the results proved in [13] and [20] have al- 
ready been applied in a variety of frameworks: see, for example, [9] for applications to 
self- intersection local times of fractional Brownian motion and [12] for an application to 
high-frequency CLTs on commutative groups. 

As an illustration of our techniques, in Section 4 we focus on generalized Volterra 
processes associated with independently scattered random measures. These processes 
have the form 

¥1"= / uh{t,s)N{du,ds), t>0, (5) 
Jr Jr 

where A'^ is a suitable Poisson measure over M x K and h is deterministic. The function 
h is called the kernel of Y'^ . When h has the form h{t, s) = v{t — s), where v has support 
in M+ ; then F'' is called a (stationary) moving average Levy process. To keep the length 
of this paper within bounds, wc only treat the specific example of a Ornstcin-Uhlenbeck 
Levy process, which is a moving average process corresponding to the case h{t,s) = 
•\/2Aexp(— A(i — s))lt>s, A > 0. These processes appear, for example, in Bayesian survival 
analysis (see James [10]), network modeling (see Cohen and Taqqu [4] or Wolpert and 
Taqqu [26]) and finance (see Bandorff- Nielsen and Shepard [3]). We will characterize the 
asymptotic behavior of the functionals 

as T — > oo. We do so by expressing the centered random variable Vr — T^^E Jq (Y/^)^ dt 

as the sum of a single and a double integral with respect to N , and we use our results 
to prove the existence of positive constants C(T) such that 



C{T) 



1 '•^ 



'™AA(0,1), (6) 



where A/'(0, 1) is a standard Gaussian random variable. 

This approach can be used beyond the framework of Ornstein-Uhlenbeck processes. 
Indeed, further applications of our main results have been developed in the two papers 
[5] and [15], where several CLTs are proved, involving linear and quadratic functionals 
of random hazard rates in nonparametric Bayesian survival analysis. Random hazard 
rates can be thought of as the mathematical representation of the 'instantaneous risk' 
associated with the average lifetime of a given population. In a Bayesian framework, they 
arc modeled in terms of moving averages of the type (5). In particular, the results proved 
in [5] and [15] can be used in prior specification procedures. See Section 4.2 for more 
details on this subject. 

As discussed below and in [5] and [15], one of the advantages of our techniques is that 
the conditions one has to check in order to prove that (6) holds can be expressed directly 
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in terms of the kernel h appearing in (5) and do not make any explicit use of (asymptotic) 
dependence properties of the process Y^^ (such as mixing). 

Observe that the theory developed in this paper also applies to a wider class of 
quadratic forms associated with Y^^ such as 

Y^^Yt''a{ds,dt) and ^ m.^^Y^Y^'', 

l<j;Vi<n 

where a and m arc, respectively, a measure and a real matrix. This kind of application, 
which will be the object of a separate study, should be compared with the findings 
contained in [1], [2] and [6]. 

The paper is organized as follows. In Section 2, we define multiple stochastic integrals 
and discuss some of their basic properties. In Section 3, we focus on sequences of single and 
double Poisson integrals, for which a general central limit theorem is proved. In Section 4, 
we describe some applications involving Ornstein-Uhlenbeck processes, random hazard 
rates and Bayesian survival analysis. 




2. Definitions and preliminary results 

Throughout the paper, (Z, Z, fi) is a Borel measure space, with /i non-atomic, a- 
finite and positive. We define the class as = {B G Z : ii{B) < oo}. The symbol 
N = {N{B) -.BgZi^i} indicates a compensated Poisson random measure on {Z,Z) with 
control /i. This means that is a collection of random variables defined of some probabil- 
ity space (f2, JF, P), indexed by the elements of and such that: (i) for every B,C G Zf^ 
such that i? n C = 0, N{B) and N{C) are independent; (ii) for every B G Zf^, 

where ^{B) is a Poisson random variable with parameter /i(i?). Note that properties (i)~ 
(ii) imply, in particular, that N is an independently scattered (or completely random) 
measure (see [19]). For every deterministic function h e L^{Z,Z,fi) = L^{fJ.), we write 
N{h) = h{z)N{dz) to indicate the Wiener-Ito integral of h with respect to N (see 
[22]). We recall that for every h € L'^{^), N{h) has an infinitely divisible law, with Levy- 
Khinchine exponent (again see [22]) given by 

V'(/i,A)= / exp(iA/i(z) - 1 - iA/i(z))^(dz), AgM. (7) 
J z 

Also, recall the isometric relation: for every (7, /i G L'^in), E[N{g)N{h)] = h{z)g{z)ij,{dz). 

Although our main results involve exclusively single and double integrals with respect 
to TV, we will often need to deal with multiple integrals of higher orders. To this end, we 
briefly recall the definition, as well as some basic properties, of general multiple Poisson 
integrals. Fix n > 2. We denote by L^(^") the space of real- valued functions on Z" that 
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are square-integrable with respect to and we write to indicate the subspace 

of composed of symmetric functions. We denote by Sn and Sn, respectively, the 

vector space generated by simple functions with the form 

f{zi, Zn) = Isi {zi) ■ ■ ■ 1b„(z„), (8) 

where Bi, . . . , i?„ are disjoint subsets of Z, and the vector space generated by the sym- 
metrizations of the elements of Sn- If / is as in (8) and / G 5^ is its symmetrization, we 
define /«(/) as 

IM) = N{Bi)---N{Bn). (9) 

The random variable /„(/) is the multiple Wiener-Ito integral of order n, of / with 
respect to N. Now, note that fi is non-atomic so /U," does not charge diagonals. This 
implies that for every n > 2, Sn is dense in and, consequently, the domain of the 

operator /„ can be extended to by continuity, due to the isometric formula, true 

for every m, ti > 2, f & Sn and g £ Sm, 

E(/n(/)/^(.9)) = n!(/,,9)i2(^„)l(„=„). (10) 

We also use the following conventional notation: /i(/) = N{f), f E L?{^j)', Inif) = Inif), 
f G L'^{fi^'^), n> 2 (this convention extends the definition of Inif) to non-symmetric 
functions /); Io{c) — c, c e R. 

Before stating the main result of the section, we recall a useful version of the multiplica- 
tion formula for multiple Poisson integrals. To this end, we define, for q,p> 1, / G L^{ijP)^ 
g G L'^ifi'^), r — 0, . . . ,q /\ p and I = 1, . . . ,r, the (contraction) kernel on which 
reduces the number of variables in the product fg from p + q to p + q — r — I as follows: 
r variables are identified and, among these, I are integrated out. This contraction kernel 
is formally defined as follows: 

/*r .9(71) ■ • -nr-hh, ■ ■ ■,tp_r, Si, • ■ ■ , Sq-r) 
f{zi, . . . ,Z;,7i, . . . ,7r_i,tl, . . - jtp-r) 

X g{zi, . . . , z;,7i, . . . ,7r-i, Si, ... , s,_,.)/^'(dzi • • -dz;) 
and, for / = 0, 

f *r .9(717 • ■ • ,7r,il, • ■ ■,tp-r,Si, . .., Sq-,.) 

(11) 

/(7I: . . . ,7r,tl, ■ • ■ ,tp-r)g{-fl, . . . ,7r, Si, . . ■,Sq_r) 

so that g{ti, • • ■ ,ip,si, . . . , Sg) = /(ti, . . . ,ip)g(si, . . . , Sq). For example, if p = g = 2, 
/*°g(7,t,s) = /(7,%(7,s), f*lg{t,s)= f f{z,t)g{z,s)fi{dz) (12) 
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(13) 

f*29= / / fizi,Z2)g{zi,Z2)fJ.idzi)fi{dz2). 



JZ Jz 

The following product formula for two Poisson multiple integrals is proved in, for 
example, [11] and [23]: letting / S L^{^p) and g G Ls(/x'^), p,q>l, and further supposing 
that f g £ i^(^P+''^'"^') for every r = 0, . . . ,p A g and I = 1, . . . ,r, we have further 
supposing 

W)i,i9)=T.^^-(l) (^)e(;) w.-k7^), (14) 

r=0 ^ ^ ^ ^ i=0 ^ ^ 

where the tilde (~) stands for symmetrization, that is, 

/ *5- 3(2^1, ■ • ■ , Xq+p-r-l) = (q +p — r — Z)! ' ■ ' ■ ' •'^<^i<l+P-r-l))j 

where a runs over all {q + p — r — 1)1 permutations of the set {1^ . . . ,q + p — r — I}. 

In what follows, we will systematically work under the assumption that there exists a 
collection of subsets {Zt :t G [0, 1]} C Z such that Zq = 0, Zi = Z, Zg C Zt for s <t and 
the following continuity condition is satisfied: for every g £ L^{fi) and every te [0, 1], 

1™/ g{z)^i(dz)= I g{z)^i(dz). (15) 

Observe that (15) is easily satisfied when Z is a Euclidean space. For instance, when 
Z — and n is equal to the Lebesgue measure, one can take Zt = [O,log(l/1 — t)]^ 
(t e [0, 1)) and Zi = R^. Also, note that for every t S [0, 1], the operator 

^t:L\ii)^L'{p):f^Titf = nz, (16) 

defines a projection operator. Since the family of projections tt = {ttj :t S [0, 1]} is non- 
decreasing and continuous, one says that tt is a continuous resolution of the identity (see, 
e.g., [28]). This concept in central in the general theory developed in [17]. In what follows, 
we will use the following notation (introduced in [19]): for every z, z' € Z , one writes 

z'^^z (17) 

whenever there exists t G [0, 1] such that z' G Zt and z G Zt - As an example, consider the 
case Z= [0,1]^* {d>l) and Zt = [0,t]'', and fix z ^ {z^^\ . . . , z^'''>) G (0,1)'*. Then z' ^„ z 
if and only if z' G [0,z)'^, where z = max(z^^\ . . . , z^'''). In particular, in the case d=l, 
one has that z' -<7r z if and only if z' < z. 

Given a kernel / G Ll{ii^) and z G Z, wc will write /!(• z) to indicate the applica- 
tion 

Z^R:y^f{y,z)liy^^zy, (18) 
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note that /!(• z) G L^{n)- As an example, again consider the case Z = [0, l]'^ and Zt = 
[0,t]'^, t G [0,1]. Then for every fixed z = (z^^), . . . , z^''') S (0,1)'', the function /!(• z) 
is the apphcation 

[0,lf^R:y^f{y,z)liye[0,zf), 
where z = max(z(^\ . . . jZ^'')). 



3. Main results: CLTs for single and double Poisson 
integrals 

In this section, we apply the results proved in [19] to study the asymptotic behavior of 
a sequence of random variables of the type 

Fn = Wn), n>l, (19) 

where /„ G L^^iJ?). In particular, our starting point is the following result, taken from 
[19]. 

Proposition 1 (See [19]). Consider the sequence {i^„:n>l} in (19) and for every 
fixed z <E Z , define /r,l(- -<7r z) £ L^{p) according to (17) and (18). Suppose that for 
every A G R, 

/ exp(iA2iV(/„l(-^,z))-l-iA2iV(/„l(-^,z)))/^(dz) ^ - — , (20) 

J Z n^oo z 

P law 

where — + denotes convergence in probability. Then Fn — > A/'(0, 1), where Af{0, 1) denotes 
a standard Gaussian random variable. 

Note that condition (20) is quite difficult to verify, since it involves a continuum of 
non-trivial transformations of the kernels /„ (one for every z G Z). However, in what 
follows, we shall show that Proposition 1 can be used to obtain neat sufficient (and 

sometimes also necessary) conditions on the kernels /„ for the CLT Fn '-^ A/'(0, 1) to 
hold. In particular, these conditions do not involve the resolution of the identity tt and 
are exclusively expressed in terms of the kernels {fn}- Our results will apply to sequences 
of kernels satisfying the following assumption. 

Assumption N. The sequence fn, n>l, in (19) verifies: 
(N-i) (integrability condition) Vn > 1, 

1/2 

jJn{z,■f^l{dz)€L\^l) and |^ /„(z, OVCd^)! e L\^Ji)■ (21) 
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(N-ii) (normalization condition) as n—s-oo, 



2 / / fn{z,zTl^{dz)fi{dz')^l; 



(22) 



Jz Jz 



(N-iii) (fourth power condition) as n—> oo, 



f f /„(^,z')V(d^)M(d2')->0 



(23) 



Jz Jz 



(this implies, in particular, that /„ e _L^(/x^) ). 

Remarks. (1) Suppose there exists a set B, independent of n, such that IJ.{B) < oo and 
for each n, /„ = fn^B, a.c.-d/i^ (this is true, in particular, when /i is finite). Then by the 
Cauchy-Schwarz inequahty, if (23) is verified, (/„) must necessarily converge to zero in 
LKfi'^). To get more general sequences (/«), we need to suppose IJ.{Z) = oo. 

(2) As shown in [13], condition (23) is not required to obtain CLTs for sequences of 
double Wiener-Ito integrals with respect to Gaussian processes and hence, in this case, 
it is not necessary to suppose fJ,{Z) = oo. Thus, in [13], the setup of a Gaussian measure 
on [0, l]'^ {d> 1) with a finite control measure was considered. 

(3) Assumption N is satisfied, for example, by a properly normalized sequence of uni- 
formly bounded functions, with supports 'slowly converging' to Z . For instance, consider 
a sequence gn G LI^^j^) such that for n>l, |(7n(-,--)|^c<oo (c independent of n) and 
the support of Qn is contained in a set of the type i?„ x B„, where < /i(i3„) < oo and 
/i(i?„) — > oo. Then if 



the sequence /„ = /i(-B„) n > 1, verifies Assumption N. Indeed, since |/„| < 

c^(B„)"\ 



The following central limit theorem is one of the main results of the paper. 

Theorem 2. Define the sequence Fn = l2(fn) o.nd /„ G L'^{^'^), n>l, as in (19), and 
suppose Assumption N holds. Then 







/n*?/nei^(At^) and fnAfneLliy?) 
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for every n>l and, moreover, 

1. z/ 

/ri*i/«^0 inLl oifi^) and /„ +2 /« ^ in L^{^), (24) 

then 

i^,/^AA(0,l), (25) 

where Af{0, 1) is a standard Gaussian random variable; 

2. if Fn G L*{V) for every n, then a sufficient condition to have (24) is that 

E(f4)^3 = E[AA(0,1)4]; (26) 

3. if the sequence {i^^:7i> 1} is uniformly integrable, then conditions (24), (25) and 
(26) are equivalent. 

Remark, (a) As already indicated in the Introduction, point 3 can be rephrased by 
saying that, if Assumption N is satisfied and {F^ ■n>l} is uniformly integrable, then 
F„ ^-^ AA(0, 1) if and only if E(F4) ^ E[AA(0, 1)*]. This result should be compared with 
the usual 'method of moments' for sequences of random variables. Sec, for example, [24] 
and [25]. 

(b) We recall (see, e.g., [8], page 355) that {Fn'-n> 1} is uniformly integrable if and 
only if 

lim supE[i^^l(j,4>j^/)] =0. 

M-.oo„>i 

(c) While the statement of Theorem 2 does not involve any resolution of the identity 
TT, the objects appearing in (16) and (17) will play a crucial role in the proof. 

(d) Observe that 

\\fn*l fn\\h^^2^ = JJ^(^Jjnia,z)fn{b,z)fi{dz)^ ^,ida) ^,idb) , (27) 

\\InAfn\?m^)= j^l^jjn{a,zfti{da)^ ^Ji{Az) , (28) 

ll/n*?/«|li.(^3)= / / /(/„(«, z)/46,z))V(d^)Ai(da)M(d6), (29) 

(e) Let G be a Gaussian measure on [Z, Z) with non-atomic control fi and for n > 1, 
let Hn = l2{hn) be the double Wiener-Ito integral of a function /i„ G Ll{y(^). In [13] 
Theorem 1 it is proved that if the normalization relation 2||/i„||^ — > 1 holds and regardless 

of Assumptions (N-i) and (N-iii), the following three conditions are equivalent: (i) iJ„ 
A/'(0, 1); (ii) E(H,t)^3; (iii) hn*\ —>0. Also, note that [13] Theorem 1 applies to 
multiple integrals of arbitrary order. 
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(f) A sufficient condition for the uniform integrability of (F^) is clearly that 
sup„E(i^^+'^) < oo for some e > 0. Note that in the Gaussian framework of [13] Theorem 
1 the uniform integrability condition is always satisfied. Indeed, by noting iJ„ = /^(ft,„) 
(n > 1), the sequence of double integrals introduced in the previous remark, for every 
p> 2, there exists a finite constant Cp such that sup„ E(|iJ„|P) < Cp(sup„ E(i7,^))P/^ < oo, 
where the last relation follows from the normalization condition E(if^) = 2||/i„|p — > 1. 

Proof of Theorem 2. Since 

fn*l fnit:S) = / /„(s,z)/„(i, z)^(dz) 



and / £ the relation /„ ★} /„ e is a consequence of the Cauchy-Schwarz 

inequality. On the other hand, by (12), 

(/n^?/n(7,i,5))V(d7,dt,ds)= f ( f fnh,tffi{dt)x f , sf fi{ds)] fi{dj) 

Jz\Jz Jz J 

^ 2 

^2 , 



/n(7,s) /^(ds)j ^(d7)<oo, 

due to part (N-i) Assumption N, so /„ ^\ /„ G 

(Proof of point 1.) According to Proposition 1, to prove point 1, it is sufficient to show 
that, when Assumption N is verified, condition (24) implies that (20) is satisfied. To this 
end, we write 

' (exp(iA2iV(/„l(- z))) - 1 - iA2Ar(/„i(. z)))/i(dz) 



— / (2iV(/„l(--<.^)))Xdz) 



(30) 



exp(iA2iV(/„l(-^,z)))-l 



- iA2iV(/„l(. z)) + ^(2iV(/„l(. z)yf 



and we shall show that under the assumptions of Theorem 2, ?7„ ^ — ^ and Vn^O- To 
show that 

(2iV(/„l(-^.^)))V(dz)^l (31) 

and hence that [/« ^ — we apply (14) in the case p = q = l and obtain 
((2iV(/„l(- z))ff 
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(32) 



^ 4 / ,Uiz,xyi{x <^ z)^l{Ax) + 4/i(/„(z, <^ z)) + 4/2(g„(z; •, ••)), 
where g„(2; ;•,••) G is given by 

= /n(2,a)/„(z,6)l(a^^^)l(b^^2). 



(33) 



The three terms in (32) correspond, respectively, to the terms (r = 1, Z = 1), (r = 1, Z = 0) 
and (r = 0, Z = 0) in (14). We deal with each term in (32) in succession. For the first term, 
observe that, due to [19] Corollary 4 and the symmetry of /„, 



fniz,x)'^fj.{dx)n{dz) 



z JZ 



/„(z, x)^[l{x z) + l{z x)]ii{dx)^i{<\z) 



Z JZ 

2 / fn{z,xfl{x z)fi{dx)fi{dz) 
J z J z 



and therefore, thanks to Assumption N 
A2 



4 / fniz,x) l{x^^ z)n{dx)fj,{dz) 



iz JZ 

A2 - 



2 / / jr,{z,xY^l{dx)^{dz) 
Z J z 



2 ■ 



(34) 



For the second term in (32), one has, by applying Fubini's theorem and the Cauchy- 
Schwarz inequality. 



E 



|/i(/„(z,-)'l(-^.^))|M(dz) 



Z2 



E[|/i(/„(z,.)'l(-^.z))||/i(/„(z',-)'l(-^^ z'))\\li{dz)ii{dz') 



< / E[/i(/„(z, ■)■'!{■ z)ff^K[h{Uz', -fli- <^ z'))Y^V(dz')A*(d^) 

JZ2 



Ai(da)/„(z,a)''l(a z) 



1/2 



/i(dz) > < oo 
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since Part ( N-i) of Assumption N is satisfied. It follows that by once again applying 
Fubini's theorem, 



E 



E[/i(/„(z, <^ z))h{Uz', -fii- z'))]^l{dz)^i{Az') 



Z2 



(35) 



fn{z.,xfl{x z)^{dz) 



/i(dx) 



< 



fn{z,xfll{Az) 



/i(dx) — * 0, 



due to (27) and (24). Hence, the integral of the second term in (32) tends to in proba- 
bility. Now, consider the third term in (32), and observe that by a Fubini argument, and 
by (33), 



E(^y_^|/2(g„(^;-,-))lM(dz) 

= / E(|/2(5„(z;-,-))/2(5n(2';v))l)M(dz)M(dz') 

< l^jniMz--.--)fY'^is{dz)j 



fn{z, affn{z, 6)^l(a^^2)l(h^^2)/i(da)/i(d6) 



1/2 



//(dz) 



iz J z 

From this, one deduces that 



fn{z,afl(a^^^)^i{da)^i{dz)^ << 



Z2 



l2i9n{z;-,-))fJ-{dz) 

nhigniz; ; ■■))l2{9n{z'\ ; ■■)Mdz)^iidz') = 2||/l„||i.(^.) 

where, thanks to (33), /i„ G L^(/i^) is such that 

/i„(a,6)= / /„(z,a)/„(2;,6)l(a^,^)l(6^„2)Ai(dz). 
Jz 

We now want to show that /„ -k\ /„(a, b) = fn{z, a)fn{z, b)^{dz) implies that ft,„ — > 
0. To do this, we start by observing that, a.e.-^^(da, d6) and thanks to Corollary 4 in 



(36) 
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[19], 

As a consequence, by noting (for fixed z) 

(z a V 5) = [{z -<7r a) n (z b)] U (a -(^ z ^„ h) [J [h z a), 

(a V -<7r z) = (a -<7r z) n (6 -<7r z), 

(a V6-(^ zV z') = [(a V6-<^ z') n (z z')] U [(a V6-<^ z) n (z' z)], 
we obtain that 



Z2 



fn{z,a)l(^a-<^^-)u{z-<^a}'i-{z-<^b)u{b-<^z)fn{z,b)fi{dz)^ M^(da,d6) 
fn{z, a)/„(z, 6)(l(^^^avb) + l(av6^,2))M(dz)^ fi^{da, db) 
/„(z,a)/„(z,6)l(2^^av6)M(dz) ) M^(da,d6) 



\Jz 



2 

2/ 



/n(2,a)/n(2,&)l(aVh^,z)M(dz)j jj. (da, d6) 
+ 2 J ^(^J f7i{z,a)fn{z,b)li^avb-<,z)lJ'{(iz)^ 

X /n(^',a)/n(2;',^)l(z^„avfe)M(dz')^/i^(da,d6). 

We now note 

(a z A z') = (a -<7r z') n (a z), 
so that, by a Fubini argument, 

/n(2:,a)/n(2,6)l(avh^,^)M(dz) ) /i^(da,d6) 



Z2 VJz 



Z2 \Jz 



/n(^:, a)fn{z', a)l^a^^zAz')fJ'{da)^ M^(dz, dz') 



and also 



2 J ^(^j fn{z,a)fn{z,b)l^aVb^^z)K'^z)^ 
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X (y^j /n(z',a)/„(z',6)l(;,/^^avb)/^(d2;')^M^(da,d6) 
fn{z,a)fn{z,h)l 
X (^J /n(^;',a)/n(2',&)l(z'A^^,av(,)M(dz')^M^(da,d6), 



\Jz 



so that the relation 



l(aV6-;^2Vz') l(2'Az^„aV6) — l(z' Az^„a,b^„zVz') + l(z'Az^„a-<„zVz') l(6-<^zAz') 

+ -'-(z'Az^„6^„zVz') l(a^„zAz') 



gives 



Z2 \Jz 



/„(z,a)/„(z,6)l(av&^„z)M(d2;) 
X (^^/„(z',a)/„(z',6)l(, f,)/x(d2;')^M^(da,d6) 

= y"^ (^y" /n(2:,a)/n(z',a)l(^Az'^„a^,zVz')Ai(da)^ pt^(d2,d2') (41) 

+ 2 J ^(^J fn{z,a)fn{z',a)l(^^^,^^a^^^y^,)fi{da)^ (42) 

X /n(z,a)/„(z',a)l(Q^^^Az')M(da)^Ai^(dz,dz'). 

After making the change of variables (z, a, 6) — > (a, z, z') in (38), observe that the terms 
(38), (40), (41) and (42) are integrals of terms of the form {A + Bf , A^, B'^ and 2AB, 
respectively, whose sum therefore equals 2{A + B^, yielding 

(/n*}/n(a,fo))'/i'(da,d6) 

(43) 

/„(z,a)/„(z,6)l(3^^avb)M(dz) ) ^^(da,d6). 



'Z2 \J z 

Since /i„ (as defined in (36)) is such that 



Z2 



hn{a,hf^i^(Aa,db)^j^(^j /„(z, a)/„(z, 6)l(av6^„z)M(dz)^ /i^(da,d6) 



and 



Z2 
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1 ^ 

f7i{z,a)fn{z,b){l(^avb^^z) + ^iz^^ayb)) K^z) //^(da,d6), 

relation (43) gives the implication: if /„ -kl /„ — > in then 

K^O inL^ifi^). (44) 

This last result, combined with (34) and (35), implies that the sequence C/„, n > 1, as 
defined in (30), converges to — ^ in probability. 

To show that Vn ^ 0, observe that | exp(iA.T) — 1 — iA.T + |A^x^| < |Axp/6 and, conse- 
quently, by Cauchy-Schwarz, 

|K|<^ / |2iV(/„l(-^. ^))lV(dz) 

1/2 1/2 ^^^^ 

< ^ (^1^ |27V(/„1(- z))\\{dz)^ iyj^ |2iV(/„l(. z))| V(dz)) 



Since the first part of the proof implies that under (24), (/^ |2A^(/„1(- z))p^(dz))^/^ — > 
1 , to conclude the proof of point 1 it is sufficient to show that under Assumption N and 
(24), 

|2iV(/„l(.^^z))|V(dz)^0 

in L^(P). To do this, one can use (32) and the orthogonality of multiple integrals of 
different orders to obtain that for any fixed z, 

E[(2iV(/„l(. <^ z))f] 

= E[((27V(/„l(--<. z))ff] 



= 16 / fn{z,xYl{x<^z)^l{Ax)\ +16/ U{z,■fl{■<^z)^Ji{da) 



iz 

+ 32 / fn{z,affr^{z,bfl(a^^^)l^^h-<^^)^?{da,db) 



and therefore 



E / \2N[fr.^■^,z))t^l[Az) 



= / E|2iV(/„l(-^,z))|V(d2) 



<16 j^[^jjn{z.xfti{dx)] M(dz) 
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+ 16 f / /„(z,a)V(da)Md2) + 32 f ( f U{z , xf ^i{dx)) /i(dz) 
JZJZ Jz\Jz J 

since Assumption N and (24) arc in order. This concludes the proof of part 1. 

{Proof of point 2.) To proof point 2, we use the product formula expansion (14) (from 
the term with r = to the terms with r = 2) to write 

Fl = Wnf = WZ^fn) + ^WT^fn) + ^h{fn A In) 

+ 2hifn 4 fn) + 2/l(/„ A fn) + 2|| i.(^2) 

and observe that since Assumption N holds and fn*^fn{a,b) ^ fn{a,b)'^ (by (ll))j 
hifn *2 /") ^ L^(P) by (23) and therefore the assumption E(F^) — > 3 implies that 

E[iF^-2h{fn4fn)f] (46) 

= E[(/4(/r^„) + AhifT^fn) + 4/2(/„ /„) + 2/i(/„ *i /„) + 2||/„||i.(^,.))'] 

^ 3. (47) 
Now, due to (46), 

E[(F2-2/2(/„*° /„))'] (48) 

= EihifT^f^f) + m&Hf^Inf) + E(16/2(/„ fnf) (49) 
+ E(4/i(/„*i/„)') + (2||/„||ij(^.))' 



There are no cross terms because the multiple integrals have different orders and 
hence are orthogonal. The most complicated term in the square bracket in (49) is 

E(/4(/„ ★(] /„)^). Since we are dealing with second order moments, the computations 
are as in the Gaussian case. We therefore obtain, using, for example, [20], formula (2), 
page 250, that 

nhif^fnf) = 4!||/^n|li2(^.) - 2(2||/„||i.(^,.))' + 16||/„*} /„||i.(^.). 

As a consequence, (48) equals 

[2(2|l/nllil„(^^))' + + 16 X 3!|l/„*? /„|li2(,,3) 

+ 16 X 2||/„*J +4|l/„*i /„||i.(^) + (2||/„||i.(^.))' (50) 

Since (50) converges to 3 by (47) and 2||/„|||2(^2) — > 1 by Assumption (N-ii), we conclude 
that ||/„*i /„||i2(^2) and ||/ti*2 /nlli2(^) ^ 0, thus proving point 2. 
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{Proof of point 3.) If F„ A/'(0, 1) and (F^) is uniformly integrable, then necessarily 
¥.{F^) E(A/'(0, 1)*) = 3. We had (26) ^ (25) ^ (24) and we just showed that uniform 
integrability and (24) imply (26), proving the equivalence of the three statements under 
uniform integrability. □ 

Example. We exhibit an elementary example of a sequence /„ € Ll{iJ?'), n > 1, verifying 
conditions (21), (22), (23) and (24). As discussed below, this example involves a sequence 
of i.i.d. random variables of the type N{B). As such, it could be alternatively worked 
out by means of a standard application of the usual CLT for i.i.d. random variables with 
second moments. However, it provides a first and simple illustration of our techniques, 
also showing that our results are consistent with the classic limit results of probability 
theory. More sophisticated examples are discussed in Section 4. Let Bj, j >l, be a 
sequence of disjoint subsets of Z such that fJ.{Bj) = 1, j > 1, and set 

Bl^={ix,y)eB,xB,:x^y}, j>l. 

Note that since fi is non-atomic, ^"^{Bq j) = iJ.^{Bj x Bj) = 1. For n > 1 and {x,y) E Z^, 
we define /„(x, y) = (2n)^^/^ i^^ v)- Of course, /„ e Ll{pL^), by definition, and 

we shall prove that {f„) also satisfies (21), (22), (23) and (24). Indeed, fn{z,-)'^fJ.{dz) 
= 2n-^j:]l,lB,i-)eL^iti) and 

n 

2||/„f = 2(2n)-i^/i'(So',) = l, 
so (/n) satisfies (21) and (22). On the other hand, 

and therefore (23) is satisfied. Finally, 



A*'(da;,dy) = -^^0 
4n 



and 



\ 1/2 

/„(z,x)V(d2:) j ^(da;) < V^/2 < oo, (51) 
fn{z,xfti{dz)] ^l{dx) = l/An^O (52) 



2 

2/ 



f„ix,z)fn{y,z)fi{dz) ] fi {dx,dy) 
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thus yielding that (/„) satisfies (24), by (27) and (28). To see (51), recall that fj.{Bj) = 1 
and write 

(^1^ Uz, x)V(d^)) ^ fi{dx) < (^(j B,^ 1^ 1b, (x)^ Kdx) 



3/2 



(we used the fact that (X]j=i (2;))^^^ = {x))- As anticipated, since (due to, 

e.g., (14)) 

71 

Wn) = 2-i/2(iV(i?,)2 - N{B,) 1), 

i=i 

the central limit result hifn) '-^ A/'(0, 1) can be verified directly by using a standard 
version of the central limit theorem, as well as the fact that N is independently scattered 
and the -Bj's are disjoint, with fJ.{Bj) = 1. 

By combining the results of [19] and Theorem 2, we can also characterize the joint 
(weak) convergence of a single and a double Poisson integral toward a bivariate Gaussian 
vector. 

Theorem 3. (A) Consider a sequence 

Gn^hign)^ [ gn{z)N{dz), 71>1, (53) 
Jz 

where gn G L'^iP') H L^{^), and suppose that as n 00, 

Ilff«lli2(^)^l and ^|g„(z)|V(dz)-.0. (54) 

Then Gn —5' X , where X is a centered standard Gaussian random variable. 

(B) Consider a sequence Fn — hifn), 'n> 1, with /„ G L^q(/x^) as in (19), and a 
sequence Gn = Iiign), 'n>l, as in part (A) above. Suppose, moreover, that: 

(i) the sequence (/„) verifies Assumption N and satisfies condition (24); 

(ii) the sequence (gn) satisfies (54). 



Then as n — s- 00 , 



(K,G„)'™(^,^'), (55) 



where X, X' are two independent, centered standard Gaussian random variables. 
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Remark. Roughly speaking, Theorem 3 tells us that for sequences of random variables 
such as {Fn) and (Gn), the (componentwise) weak convergence of F„ and G„ toward 
a Gaussian distribution always implies the joint convergence of the vector (i^„,G„). As 
proved in [20], an analogous property holds for vectors of multiple Wiener-Ito integrals 
with respect to general Gaussian processes. 

Proof of Theorem 3. (A) By using (7), for every A S M, 
E[cxp(iAG„)] 



; exp 



exp 



[exp(iAg„(z)) - 1 - iXgn{z)]^l{dz) 



(^-yllSrJIi^f^)^ cxp 



exp(iAg„(z)) - 1 - iXgn{z) + —gnizf 



/i(dz) 



To conclude, observe that by (54), ||,9n Ili2(^-) 1 and the integral in the second expo- 
nential is bounded by 



cxp(iAg„(z)) - 1 - \\gn{z) + —gn{zf 



|g„(z)|V(dz)~>0, 



thus yielding E[exp(iAG„)] ^exp(— A^/2), as required. 

(B) To prove (55), it is sufficient to show that for every a,(3 & 



(56) 



We shall prove relation (56) by using [19] Theorem 7, as well as some estimates appearing 
in the proof of Theorem 2. To do this, define X^(K) = h{z)N{dz) = Ii{h) for every 
h G L'^{Z, Z, ^) = Sjf^ and consider the random field 

X^^{X^{h):heSi,}. 



Note that X^ belongs to the class of random fields studied in [19]. Now, for every n, 
denote by tt the resolution of the identity over the Hilbert space given by 

Trth{z)^lz,{z)h{z), h e Sjf„t G [0,1], 

where the sets Zt appear in (15), and further define 

ui^'^\z) = 2a/i(/„(.,z)l(. z)) + /?5„(z) 

ah^{f,,){z) + f3gn{z) e L^(f)^,X^), 

where the class L^{Sjfi,X^) of 7r-adapted random functions is defined in Section 3.2 of [17] 
and the process z i-^ K{fn)[z) G L^(^^,X-) is defined as K{fn){z) = 2/i(/„(-, z)l(- 
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z)). By using [19] Proposition 5, we can write 

N N N 

where the generahzcd stochastic integral is defined in Section 3.2 of [17] and the 

JV 

second equahty is a consequence of the Unearity of the operator Jx^- From [17] Theorem 

N 

7, in we deduce that (56) is proved, once it is shown that 



(57) 



To prove (57), we adopt the same strategy as in the proof of Theorem 2. that is. we 
write 

" [exp(iAu("-'5)(z)) - 1 - iA«(-«(z)]/i(dz) 



{u\-^P\z)f^,{dz) 



/i(dz) 



and show that under Assumptions (i)-(iii) in the statement, An —^{a^ + /3^) and 
Bn 0. We start by considering i3„, writing 



\BJ''' 



cxp 



(iAu("'«(z)) - 1 - iA^("'«(z) + -{u^^^'\z)f 



^l{Az)\ 



1/3 



< 



61/3 



|z.(r'«(z)|V(dz) 



< 



Ml 

6 



1/3 



1 /3 

^ (/^ |/i.(/„)(^)l V(d;.)) + ^ (/^ l<?n(^)l V(dz) 



61/3 



|a/i^(/«)(^) + /?g„(z)|V(dz) 

1/3 



1/3 



Now recall that in the proof of Theorem 2, we already verified in (45) that under As- 
sumption N and (24), 

' \K{fn){z)?^l{Az)^0. 

Since (/^ \9n{z)\^ ^J^{'A^z)Y/^ (due to (54)), we deduce that B„ ^ 0. To prove the 
convergence of the sequence A„, we write 

(^i"-«(.))V(dz) 
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+ 2ap / /i,(/„)(z)g„(z)M(d^). 



In the proof of Theorem 2. we showed that when Assumption N and (24) hold, 
/i^(/„)(z)^ X n{dz) 1 (see (31)). Since g„(z)^^(dz) ^ 1 by (54), to prove that 
An ^ — 4-(q:^ + /3^), wc have only to show that 



K{fn){z)gn{z)^Ji{(lz) 0. 

To this end, we use the definition of h.^{fn) and apply Fubini's theorem to obtain that 



E 



\K{fn)iz)gn{z)\^J,{dz) 



< 



fn{a,zyi{a-(^ z)fi{da) / f„{b, z'yi{b z')i^{db) 

Z2 Uz 

X \9n{z')gn{z)\fl{dz)fl{dz') 



1/2 



.1/2 y 

I Kdz)j 



|.g„(z)| / /„(6, 2)^1(6 ^,z)A.(dfe) ^iidz) 



< / .9„(z)V(dz)x / / Ub,z')H{b^^z')fi{db)\fi{dz')<(x,. 



Z IJZ 



As a consequence, by once again using the isometric properties of the random measure 
N, Fubini's theorem and Cauchy-Schwarz, 



Z2 



h^ifn)iz)gniz)fi{dz) 
E[/i(/„(., z)l(. z))5„(z)/i (/„(., z')l(- z'))gr.{z')]^l{dz)^l{dz') 
fn{a,z)l{a^.^ z)g„{z)fi{dz)) ^(da) 
/,i(a,z)l(a z)gn{z)^i{dz) 



z \J z 

X fn{a,z')'^{a z')gn{z')^Ji{dz')^^l{da) 
gn{z)gn{z') (^J fn{a, z)fn{a, z')l{a z A z')^ /"^(dz, dz') 
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X (^J ^(^J fn{a,z)fnia,z')l{a^„ z Az')fi{da)^ ^'^{dz,dz')J . 
Note that (/^^ 5«(^)^gn(2') V^(d2: dz'))!/^ = ,g„(z)2^(dz) 1, by (54). Moreover, 
fn{a, z)fn{a, z')l{a -(^ z A z')^(da)^ /x^(dz,dz') 

(^J fn{a,z)fn{b,z)l{aV b^^ z)fj.{dz)^ /i^(da,d6) 
hn{a,b)^fi^{da,db)-^0, 

Z2 



Z2 \Jz 



where we again appUcd Fubini's theorem, /i„ is given by (36) and we used (44). It foUows 

that An — ^(a^ + (3'^) and therefore that (56) holds, thus concluding the proof of 
Theorem 3. □ 



4. Applications 

This section contains two illustrations of our techniques. Both involve generalized Voltcrra 
processes, that is, random processes having the form (5), where iV is a Poisson measure 
over R X R (e.g., with control measure i'{du) ds) and ft, is a deterministic bivariatc kernel. 
In Section 4.1, we prove CLTs involving linear and quadratic functionals of Ornstein- 
Uhlcnbcck Levy processes, obtained from (5) by setting 

h(t,s) = V2Xe^p{-X{t- s))lt>s, A>0. 

See [3], [10] and [26], respectively, for applications of Ornstein-Uhlenbeck Levy processes 
to finance, survival analysis and network modeling. Section 4.2 contains a concise de- 
scription of the applications of the theory of this paper developed in [5] and [15], where 
Theorem 2 and Theorem 3 arc applied in order to obtain CLTs involving random hazard 
rates in Bayesian survival analysis. 

4.1. Ornstein-Uhlenbeck Levy processes 

Fix A > 0. We consider the Ornstein-Uhlenbeck Levy process given by 

Y/'=^V2x[ [ uexp{-X{t-x))N{du,dx), t > 0, (58) 
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where iV is a centered Poisson measure over R x K with control measure given by 
i^(dit) da;, where i'(-) is positive and normahzed in such a way that JjjU^z^(dit) = 1. We 
also assume that \u\^i'{du) < oo. Note that Y^.^ is a stationary moving average Levy 
process. We shall use part A of Theorem 3 to prove the following result involving linear 
functional of F . 



Theorem 4. As T ^ oo, 



VtJo 



(59) 



where N{0,a^{X)) stands for a centered Gaussian random variable with variance cr^(A) 
Proof. First, put the integral in the form (53). Applying Fubini's theorem, we have 



Yt^dt 



1/2 i-T 



cxp(— A(t — x)) dt 



N{du,dx). 



We need to verify (54). First 

i-T \ 2 



E 



7f i„ 



2A 



1/2 



— ) A-icxp(Ax)(l-exp(-AT)) 



i^(du) dx 



2X 





1 /.O 



2A 



1/2 



T 



— ) A ^ exp(Aa;)(exp(— A.t) — exp(— AT)) 



[exp(Ax)(l - exp(-Ar))]Mx 



^{du) dx 



2A 



-1 rT 



T 



[exp(Aa;)(exp(— A.t) — exp(— AT))]^ dx 



-^2A"\ asT^oo. 
According to Theorem 3, to prove (59) it is now sufficient to show that as T — > oo, 



1/2 i-T 



cxp(— A(i — x)) dt 



^{du) dx — > 0. 
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But, 



1/2 



exp(— A(i — x)) dt 



\ufiy{du) X 



/2A 



a:VO 
-1x3/2 



i^{du) dx 



V T 



[cxp(Ax)(l - cxp(-AT))]''da; 



(60) 



oo 
T 



+ / [cxp(Aa;)(cxp(-Ax) -cxp(-Ar))]''dx 



< / |u|^i/(du) X 
Jr 

thus concluding the proof. 



/2A- 



-1\ 3/2 



3A 



T 



•0, 



(61) 
□ 



A CLT analogous to (59) also holds in the more general case of an extended Ornstein- 
Uhlenbeck Levy process Y^'^ ^ defined as 



4>,\ 



texp(-A(t - x))(l){t,x)N{du,dx), t > 0, 



(62) 



where A > and the kernel (j){t,x), from R_|_ x M to M+, is such that < e < 4>{t,x) < 
7] < oo for some strictly positive finite constants e and rj. Indeed, in this case, one can 
easily show that the variance of the random variable Jg^y/'^dt, denoted V^^x(T), is an 
increasing function (whose explicit expression depends on the choice of A and (f>) such 
that 



a(£,r;,A)<-5^^<6(£,7y,A), 



r>o, 



(63) 



where a(e, rj, A) and 6(e, rj, A) are two positive constants not depending on T. Finally, the 
combination of the estimate (63) and of the arguments displayed in formulae (60) and 
(61) yields the CLT 



y/'^dt'™AA(o,i), 



T oo. 



We now turn to quadratic functionals. We shall now suppose that the measure is 
such that /w^V(du) < oo for j = 2,4,6, and / w^z^(du) = 1. Plainly, these assumptions 
yield that 



[m1-2A 



Theorem 5. For every A > 0, define 



iV2^(*-^V(du)da;==l. 



h^{u,x;u',x') = 2Xuu' exp{~X{t — x) — X{t — x'))l(^_oo,t]^ix,x'). 



(64) 
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Then as T — > oo ; 
1. 



dt 



(65) 



(66) 



where I2 is a double Poisson integral with respect to N , and Af and Af' are two 
independent centered Gaussian random variables with variances given, respectively, 
by and c1 = J u'^i^{du); 



2. 



where M {0,1) indicates a centered standard Gaussian random variable. 
Proof. We introduce the notation 



1 



(67) 



Jf,,T(«,x;u',x') = (u X «')^^^2^^4!^^{e^(^+^'ni -e-^^)l(.v.'<o) 



T 



-2XT 



HIt{u,x) = U 



2 l(-oo,T](a;) 



T 



{e^-'^(l 



-2T 



)l(,<o)+e2^-(e 



~2Xx 



)l(xV3;'>0)}7 (68) 
(x>0)}- 



A standard interchange of deterministic and stochastic integration yields that 



2Au2e-2A(*-^)iV(du,da:) 



dt 



= {l2{VTHx,T),Ii{VfHx,T)}^{K2{T),Ki{T)}, 

where I2 and Ii denote, respectively, a double and a single Wiener-Ito integral with 
respect to TV (observe that, since in (65) the parameter t is integrated with respect to a 
finite measure, the interchange of deterministic and stochastic integrals can be justified 
by means of a standard stochastic Fubini theorem - one can, e.g., mimic the proof 
of [16] Lemma 13, by first approximating the kernels and M2g-2A(t-x) -^y. j^gg^ns of 
piecewise constant integrands and then by using the isometric properties of single and 
double Wiener-Ito integrals). Since for every T, Ki{T) is a single integral and K2{T) is 
a double integral, the joint convergence of the vector {K2{T), Ki{T)} can be studied by 
means of Theorem 3. By using the same kind of calculations as in the proof of Theorem 
4, one easily verifies that 



l[Ki{Tf]^T [ Hlj,{u,xfv{du)dx^ 
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and 



We therefore deduce from part A of Theorem 3 that Ki(T) ^ N{0,cl). In view of part 
B of Theorem 3, the CLT (66) holds if the kernel 



Jt{u, x; u', x') = VTHx t{u, x; u' , x') (69) 

verifies Assumption N (i.e., relations (21)-(23)) and (24). Start by observing that (N-ii) = 
(22) holds (up to a different normalization) because 

2||>^T||L2((ExR)2,(di^xda;)2) 

= 2T||i7A,T||i2((RxR)2,(dl/Xd2;)2) 

2 



/ JKxR 

X e^(-+-')(e-2^(^^^') -e-2^^)l(,v.'>o)}'dxda;' 



It is therefore sufficient to verify that kernel Jt defined in (69) verifies conditions (N-i) = 
(21), (N-i) = (23) and the contraction condition (24), namely that 

JT{-]u,xfv{du)dxeL^{M.x%dvdx) Vr>0, (70) 

R 

. 1/2 

JT{-]u,xfv{du)dx\ e ^^(M X R,d!^d.T) Vr>0, (71) 



JT{u,x;u\x)'^v{du)dxv{du')dx T^oo, (72) 

II Jt *2 -^T||i2(RxR,di/d2:) ~* 0, T-^oo, (73) 

II Jt*} JT||i2((RxR)2 ^ 0, T^oo. (74) 

In view of (69), (70) and (71) are implied by the definition of H\ t- Relation (72) can 
be deduced from the relation 

T'x||^lTlli2(R,R,d.d.)-^-0. (75) 

Relations (73) and (74) are a consequence of the two asymptotic relations 

2 



X / Hx,T{u,x;u',x'yv{<iu)dx] v{du')dx- 
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(76) 



T I / Hx^T{v,z','U',x)H\^T{v,z;u',x')i'{dv)dz\ i'{du) dxv{du') dx' 

J(RxR)2 \Jmxm ' ' / 

(77) 

Note that (75)-(77) can be easily checked by resorting to the exphcit definition of Hx^t, 
as given in (68). The conclusion of point 1 now follows. Point 2 can be deduced from the 
relation 

Vr|l^^(r,^)2dt-r 

= Vf(^^£l2{h^)dt+^£ I' l^2Xu^e-^^^'~^m{du,dx) dtj, 

which is a consequence of the multiplication formula (14) (in the case p = q = 1), applied 
for every t to the variable (Y^'^)^. □ 

Remark. In the Poisson case considered here, the double and single integrals in (66) 
both converge to a non-degenerate Gaussian distribution. This situation is different when 
the random measure is Gaussian. In that case, the single integral - which corresponds to 
the contribution of the diagonal - is deterministic. 

We can also prove the following CLT for the sample variance of Y/". 

Corollary 6. With the same notation as in Theorem 5, for every A > 0, 

Proof. We simply write ^ {Y,^ - ^ Y^^ du)^ dt = ^ (Y,^)^ dt - (i/„^y,Mt)2 
and observe that by (59), Vt{^ Y/" dt^ = OviT-^'"^)- □ 



4.2. Functionals of random hazard rates in nonparametric 
Bayesian survival analysis 

We shall now discuss some applications to random hazard rate models in nonparametric 
Bayesian survival analysis. These random hazard rates are often represented as gener- 
alized Volterra processes of the kind described above. In [5] and [15], the linear and 
quadratic functionals associated with random hazard rates in some popular Bayesian 
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models are studied by means of the techniques developed in this paper. In what follows, 
we present the main elements of a Bayesian random hazard rate model, as well as some 
crucial examples and motivations taken from [5] and [15]. The reader is referred to [5], 
[10], [15] and the references therein for more details about Bayesian models of survival 
analysis. 

Let {Z,Z) be a measurable space and let /i be a cr- finite measure over {Z,Z). As 
before, we use the notation = {B G Z : ^{B) < oo}. A collection of random variables 
N = {N{B) : B G Z^} is called a non-compensated Poisson measure with control measure 
fi if there exists a compensated Poisson measure N = {N{B) .B £ Z^} (see Section 2) 
such that 

N{B) N{B) + n{B) VBeZ^. (78) 

In other words, = {N{B) : B G Z^} is a non-compensated Poisson measure if and only 
if (i) for every B,C £ Z^ such that B O C = , N{B) and N{C) arc independent, and 
(ii) for every B £ Z^, 

N{B) '=^3(5), 

where ^{B) is a Poisson random variable with parameter fJ,{B). 

The basic ingredients of a Bayesian random hazard rate model are the following: 

• a random hazard rate h, which is a positive generalized Volterra process of the type 

h{t)^ / / uk{t,x)N{du,dx), t>0, (79) 

where k{t, x) >0 and A^ is a suitable non-compensated Poisson measure over R x ]R+ , 
verifying 

lim / h(t)dt^+oo, a.s.-P: 

• a random density with support in , given by 

f{t) = h{t)cxpl^- J^h{s)ds^ = h{t)exp{~H{t)}, t>0, (80) 

where 

H{t)= I h[s)ds, t>0; 

• a sequence of positive absolutely continuous exchangeable random variables U = 
{Un '■n> 1}, representing the lifetimes associated with a given population, such that, 
conditionally on the density / in (80), U is composed of i.i.d. random variables with 
common law given by /. 
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Plainly, the initial choice of the law of h and / (which, in Bayesian terms, is the prior 
specification of the model) is completely encoded by the choices of the control measure 
of A'' and of the kernel k. Note that (80) gives the following heuristic characterization of 
hit): 

h(t) dt^F(t<Ui<t + dt\Ui >tj), 

meaning that given /, the quantity h{t)dt is the probability that the lifetime Ui (or, 
for that matter, Un) falls in the interval [i,t+ dt], conditionally on the fact that ?7i is 
greater than t. 

Popular choices for k arc the following: (a) the Dykstra-Laud kernel k{t,x) = lo<a;<t; 
(b) the rectangular kernel k{t, x) = 'i-\x-t\<T} where r > is called the 'bandwith' of the 
kernel; (c) the Ornstein-Uhlenbeck kernel k{t,x) = '\/2Aexp{— A(i — a;)}lo<a;<t (note the 
difference with (58), where the kernel is indeed •\/2Aexp{— A(t — a;)}l_oo<a;<t)- 

Popular choices for the control measure of N are: 

(i) generalized Gamma controls of the type 

li{du, dx) =T{1 — (j)^^ cxp{—ju)u^^^'^lu,x>o dudx, (81) 

where a G (0, 1), 7 > and T is the usual Gamma function; 

(ii) extended Gamma controls of the kind 

fi{du, dx) — exp{~ f3{x)u)u~^lu x>o dudx, (82) 

where /3 is a strictly positive function on R; 

(iii) Beta controls of the type 

H{du,dx) = (1 -M)'='"^)-^c(a;)l„g(o,i)lx>odudx, (83) 

where c is a strictly positive function on IR. 

Note that, in general, extended Gamma controls and Beta controls are non- 
homogeneous, in the sense that they cannot be represented in the form ij,{du,dx) = 
i'(du)dx, for some tr-finite measure v on R^. 

The crucial point is that due to the very complex nature of an object such as (79), very 
little is known about the effect that different parametric choices in (81)-(83) may have 
on the distributional behavior oi h, in particular with respect to functionals of statistical 
relevance. The idea developed in [15] and, later, in [5], is that one can always represent 
h in terms of some underlying compensated Poisson measure (via relation (78)) so that 
Theorem 2 and Theorem 3 of the present paper can be applied in order to obtain CLTs 
for linear and quadratic functionals of h. The key step is, of course, to represent linear 
and quadratic functionals as a sum of a single and a double Wiener-Ito integral, a task 
that can be easily performed by using, for example, the multiplication formula (14). It 
turns out that in most cases, the constants involved in the CLTs obtained in this way 
can be expressed very neatly in terms of the different parameters composing the control 
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measure of N and the kernel k, thus giving a rough description of the overall 'shape' of 
the hazard rate h as a function of the prior specification of the model. As argued in [15] 
and [5] , these kind of results may serve as a guide in the prior analysis since they provide 
a quite direct way to incorporate prior knowledge into the specification of the law of h. 

As an example, we present some results, proved in [15] by means of Theorem 2 and 
Theorem 3, involving linear and quadratic functionals of the rectangular kernel, under 
different parametric choices of the control measure associated with TV. The first statement 
concerns linear functionals. Recall that the random variables H{T) = /J"/i(s)ds, T> 
0, have been defined in (80); the symbol A/'(0,c) denotes a centered Gaussian random 
variable with variance c. 

Theorem 7 (See [15]). (1) Let N have a homogeneous control measure of the type 
/i(du,dx) = v{du)(lx, where the a-finite measure v is such that K^'^ = uV(du) < oo, 

i = 1,2. Let h{t) be defined via (79), with the rectangular kernel k(t,x) ~ 1|j_j.|<t-, for 
some T > 0. Then 

-^[H{T) - 2Ti^(i)T] ^ AA(0,4T2i^(2)). 

(2) Let N have a non-homogeneous control measure of the type (82), with P{x) = 
1 + x^/^. Let h{t) he defined via (79), with the rectangular kernel of bandwith one given 
by k{t,x) = lit-x\<i- Then 

1 r~ /„s .^1/2] law 



[i/(T)-4T^/i-^AA(0,4). 



(3) Let N have a non-homogeneous control measure of the type (83), with c{x) ~ x^/^ 
as X ^ oo. Let h(t) be defined via (79) with the rectangular kernel k(t,x) = l|t-a;|<i. 
Then 

^ [i/(T)-2r]^AA(0,8). 



Ti/4 



The following statement involves quadratic functionals of rectangular random hazard 
rates, under a homogeneous assumption on the control measure of N . Further results, 
involving non-homogeneous random measures, are contained at the end of [15] Section 
4.2. 

Theorem 8 (See [15], Section 4.2). Let N have a homogeneous control measure of the 
type /i(dit,dx) = v{du) Ax, where the a-finite measure v is such that K^^ = J^^ u'i/(du) < 

oo, i = 1, . . . , 4. Let h(t) be defined via (79), with the rectangular kernel k{t, x) = '^\t-x\<T! 
for some t > 0. Then 



pT 



i r h{tf At - (2TAf ) + 4t2(A-(i))') 



T 







^AA(0,ci) 
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with ci = 16T2[i^(^) /4 + TKl'>K)f> + 2T[ii:i^^] 2/3 + T^[K}:f>YKl'>]. Moreover, 



'(1)7^(3) 



(2)12, 



.2 r 7^(2)12 ^-(1)1 



1 



H(r) 



where - ^t'^\K^^^ + 8T[ii:i^^]2/3] 



dt - 2x1^(2) 



law 



AA(0,C2), 



Remark. The paper [5] continues the analysis contained in [15] by estabUshing posterior 
CLTs for Unear and quadratic functionals of random hazard rates, that is, hmit theo- 
rems involving the law obtained by conditioning on an arbitrary sample of observations 
(C/i, . . . , [/„). One of the main findings in [5] is that the posterior CLTs coincide with 
the prior ones for most models commonly used in Bayesian analysis. These results are 
then compared with another asymptotic characterization of Bayesian models, known as 
consistency (see, e.g., Draghici and Ramamoorthi [7]). 
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